Magnetic dynamics driven by the spin-current generated via spin-Seebeck effect 
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We consider the spin-current driven dynamics of a magnetic nanostructure in a conductive mag- 
netic wire under a heat gradient in an open circuit, spin Seebeck effect geometry. It is shown that 
the spin-current scattering results in a spin-current torque acting on the nanostructure and leading 
to precession and displacement. The scattering leads also to a redistribution of the spin electrochem- 
ical potential along the wire resulting in a break of the polarity-reversal symmetry of the inverse 
spin Hall effect voltage with respect to the heat gradient inversion. 



Introduction.- The discovery in the 1820s by T.J. See- 
beck that due to a temperature gradient an electric volt- 
age emerges along the temperature drop, revealed the 
relationship between heat and charge currents. The re- 
versal of Seebeck's effect, i.e. the appearance of a tem- 
perature gradient upon an applied voltage, was shortly 
thereafter confirmed by J.-C. Peltier in 1834. In addition 
to other thermo-electric phenomena such as the Joule 
heating, in magnetic fields new thermo-magnetic effects 
arise: A resistive conductor with a temperature gradient 
VT placed in a magnetic field B perpedicular to VT 
develops a potential drop normal to both VT and B. 
This phenomenon is termed the Ettingshausen effect and 
its reverse is the Nernst effect. In a magnetic material 
the anomalous Nernst effect occurs (i.e. the Nernst effect 
due to the spontaneous magnetization) which was first 
observed for Ni and Ni-Cu alloy^^- . Recently, in Refs.-^^ 
measurements of the planar and the anomalous Nernst 
effect were reported for a variety of materials includ- 
ing magnetic semiconductor, ferromagnetic metals, pure 
transition metals, oxides, and chalcogenides. A quali- 
tatively new phenomena, the Spin-Seebeck effect (SSE), 
was discovered 2008 by Uchida et alJ showing that in a 
ferromagnetic material (a mm-size NigiFeig sample) and 
in an open-circuit geometry (which is also the geometry 
studied in this work, cf. Fig. [T|) a heat current results 
in a spin current, i.e. a fiow of spin angular momen- 
tum and hence a spin voltage, even if VT is parallel to 
B (where the Nernst-Ettingshausen effect does not con- 
tribute). The spin voltage is reflected by a charge voltage 
that emerges (due the inverse spin Hall effect (ISHE)) in 
a Pt strip deposited on the sample perpendicular to VT 
(cf. Figll]). Further experiments on resistive conductors 
(— for NigiFcig), insulating ferrimagnets (LaY2Fe50i2 
inS), and for ferromagnetic semiconductors (GaMnAsiS) 
underline the generality of SSE. These fascinating effects 
are not only of fundamental importance; thermo-electric 
elements are already indispensable for temperature sens- 
ing and control and for current-heat conversion. SSE 
opens the way for thermo-electric spintronic devices with 
qualitatively new tools for energy-consumption reduc- 
tion. It is highly desirable to explore whether SSE can be 
utilized to steer localized magnetic textures, as problem 
addressed here. Theoretically, the reciprocity between 



the dynamics in the magnetic order and the heat gradi- 
ents is governed by the Onsager relations. The Onsager 
matrix were discussed from a general point of view in 
Ref.JA with a focus on the transport of charge, magneti- 
zation, and heat. In RefSi ^^d^ a thermo-magnetic meso- 
scopic circuit theory was put forward. Refi^ pointed 
out the occurrence of a thermally excited spin current 
in resistive conductors with an embedded ferromagnetic 
nanoclusters. Other recent works^^'^"^ addressed the ther- 
mally induced spin-transfer torque in a spin valves struc- 
tures whereas the phenomenological studj»i^ is focused on 
the spin-transfer torques in quasi one-dimensional mag- 
netic domain walls (DWs) by introducing a viscous term 
into the Landau-Lifshitz-Gilbert equation (LLG). 
Spin current- The microscopic mechanism for the appear- 
ance of the spin current in SSE is not yet completely 
understood, it is however an experimental fact that in 
the geometry of FiglT] the thermal gradient VT gener- 
ates a steady state spin current without a charge 
current!"— . The amplitude of is found to be deter- 
mined by2r— 
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where k is a temperature-independent transport coeffi- 
cient whose properties are discussed in^^lO; no charge 
current is generated. The purpose of this work is to in- 
spect the dynamics triggered by (eq.((T])) for the case 
where localized magnetic texture Mi^ is present in the 
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FIG. 1. A localized magnetic structure M in a ferromagnetic 
conductor of length 2L subject to the constant temperature 
gradient VT. B is a saturation magnetic field along VT. 
The pure spin current JJ is signaled by the inverse spin Hall 
voltage ^^jsjjjt; measured in a Pt strip, x and z axes are 
indicted. 
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ferromagnetic (FM) conductor (cf. Fig[T]), a problem of 
great importance and has not been addressed so far. As 
show below, the quantum-mechanical scattering of 
from M acts in effect with a spin-current torque on M 
which results in an oscillatory and a displacement motion 
of M. Upon scattering also changes. This leads to a 
redistribution of the spin electrochemical potential which 
can be measured via ISHE. 

The system under consideration is illustrated in Fig. 
[TJ Two thermal reservoirs with different temperatures T 
and T + AT create along the x axis in a FM conductive 
wire of length 2L a steady T-gradient VT and hence a 
steady-state spin current J^. This means, without know- 
ing the detail of the operators associated with SSE, these 
project the system onto a chargeless eigenstate "ipix) of 
the spin current operator j*(/c). Generally, such a state 
can be written asi^ 
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The X expectation value of the charge current 'f{k) — 
§^[{y^l)\x))^l;{x)-ipHx)Vip{x)] vanishes, i.e. = 

(here m is the effective mass). In contrast, for the spin 
current j^(fc) - ^[{Vij\x))a^i;{x) - tl^^ {x)a^Vij{x)] 
we infer 



^■^(^) = l^(cos2(/>-cos20). 
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In general, the thermal transport is ballistic^° but with 
diffusive spins, i.e., upon creating ([2]) the spin coherence 
is lifted by scattering events that randomize 4> and 9 
within [0, 27r]. Hence, the expectation value of the spin 
current vanishes on the scale of the spin-flip diffusion 
Iengtb2i.i22, i.e., r^{k) = § j'^dedcf) / {2n)'^ = 0. However, 
when the wire is magnetically polarized and driven to 
saturation by the magnetic field Bl"— , we find (ax) ^ 
but {(7y) — 0. Eq.(l2]) reads then for an exchange-split 
conductor 
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Here 6* G [0, 27r] still appears due to the residual spin 
precession and diffusion. Then we have J'^ = 0, J^(fc) = 
§ jydO/2'K = 0, whereas J|(fc) — J^ik) — hk/2m in fine 
with the experimental observation^"^*^. 

The main purpose of the present work is to investigate 
the influence of a localized magnetic, non-diffusive scat- 
terer M (where x = is taken as its central position (see 
FigllJ). M has a uniaxial anisotropy along an axis chosen 
to be z. If M(a;) has an internal structure, e.g. a non- 
collinearity, which varies on a scale larger than X = 2n/k 
(the variation scale of ([5])), one can unitary transform to 
align locally with M which introduces a weak gauge po- 
tential that can be dealt^^ with in a perturbative way us- 
ing the Green's function constructed from ([5]) (similarly 
as done i n^^i^^ ). We find M has a stronger influence if 



its range of variation is comparable to A. In this case 
the magnetic texture acts in effect as M(x) — M.oS{x), 
where the magnetic moment Mq derives from an average 
of M(a;) over its extension w. The model is realizable 
for^'^ rather then for metals. The interaction between M 
and the electron spin cr reads^. 



Hint = gM{x) ■ (T 



(6) 



where g is a local coupling constant and A/q is large 
enough to be treated classically. For Mq aligned with 
z axis, as in Fig HI we derive using eas. (|5l6p the expres- 
sion for the spinor wave function in the presence of M(a;), 
namely 
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forx<o, 

(7) 

The scattering state V's(a;) describes the spinor wave 
function in the original spin channel, which is partially 
reflected into the original-spin and the spin-flip channels, 
and also partially transmitted into theses two channel, 
which gives the complex spin reflection and transmission 
coefficients r and 
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a = gMom/kh'^. (8) 



The magnetic scattering gives rise to a short pseudo- 
circuit to the charge channels, for we find 



2a sin 9 
1 + 



(9) 



(j'l) vanishes however beyond the spin diffusion length 
after averaging over 9. Counterparts, i.e. a pure spin 
current generated by a charge current when scattered off 
a magnetic structure, are well known, e.gi^^— . The spin 
current carried by ([5]) is modified upon scattering and a 
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FIG. 2. The electric voltage Vishe generated by the inverse 
spin-Hall effect in the Pt layer as a function of the mag- 
netic scattering strength ap ~ gMom/kpf? with kp being 
the Fermi wave vector. The relaxation rate is T/Ep = 0.01. 
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none-zero J* emerges 
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for X < 0, 
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Within linear response, the spin voltage along the 
wire is fisix) = £,sx{J^{a)) where is a function 
of the elementary charge, the spin-dependent electric 
conductivity, the spin-dependent Seebeck coefficient, 
and the spin-Seebeck coefhcient of the FM wire^S. 
The quantum-mechanically averaged spin current is 
(J|(a)) — Tr{J^{a)p)^. The single electron density ma- 
trix is p = ^^l^X^k -Ef - g^-i-ir '^here Ep is the Fermi 

energy, — h^k'^ /2m + ?i^k|/2m with k|| being the 
transverse wave vector, and F is a Lorenzian relaxation 
rate due to disorder"^ Depositing a conductive strip 
with a strong spin-orbit coupling, e.g. Ft, as shown in 
FiglH ^s{x) can be imaged via the electric voltage Vishe 
generated by the inverse spin-Hall effect in Ft using the 
relation 



Vishe (a^) _ x 
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where Vishe(^) is the electric voltage measured in 
Ft in absence of the magnetic scatterer M. Explic- 
itly, VisHE(a;) = j^sx{J^{0)) with 7 being a system- 
dependent parameter—^—, determined by the spin-Hall 
angle in Ft, the spin-injection efficiency across the 
FM/Ft interface, and the length and the thickness of 
Ft wire. Due to the spin current scattering off M the 
Hall voltage Vishe losses its odd symmetry with respect 
to a reflection at x = 0, i.e. we deduce —Vishe (—2:) ^ 
Vishe (a;)- As shown in Figl5]the amount of the symmetry 
break depends on a, and can be taken in the experiment 
as an indicator of the presence of magnetic scattering 
centers. 

Magnetization dynamics.- In as much as J| is modified 
by the presence of M, the scattering triggers a dynam- 
ics of M which is usually much slower than the carrier 
scattering dynamics and can be classically treated (Mq 
is assumed large). J| acts on M with a torque that 
follows from the jump in the spin current at the point 
X = 0,T^ = J'^{0-) - J^{0+). Hence, derives from 
our quantum mechanical calculations as 



Tx — Jq 
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Both components are transversal. Ty tends to align M to 
the direction of the FM magnetization, while tries to 
rotate the moment M around the axis Bx- Equivalently, 
within our model, the spin-current torque is obtained 
from the spin density Sfj_{x) accumulated at the localized 
moment (due to interference of incoming and reflected 
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FIG. 3. Precession of the magnetic moment M (eq.((T6}) for 
different spin-current scattering strength a (cf. eq[8]) and 
Gilbert damping ag. Here Vx = dX/dt and D^/g = 5. 
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FIG. 4. Displacement of the magnetic moment (given by 
ea. (|17p ) vs. time for the same parameters of Fig|3] Vx = 
dx/dt. 



waves) as 
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[n X S{x = 0)] 
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where n is the unit vector along M, and the spin den- 
sity we obtain from Sf^{x) = ^^(x)ct^'!/'(x). Since Mq is 
assumed large (say > 5/2 fig) the spin current-induced 
magnetization dynamics can be treated with the modified 
LLG equation^i^ 



dn _ D^, 
'dt ~ 



n X e. 



9n g 
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n X S(0)] (15) 



where Dz is the anisotropy energy and ag is the Gilbert 
damping parameter—. Two motion types of M occur: 

Precession- Introducing the following magnetization 
distribution 



n= [mu (t) sin X{t),my{t),mn{t) cos X{t)] (16) 
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and propagating with the LLG equation fEg fTS]) starting 
from r7i||(0) = 1, TOy(O) = 0, and X(0) = 0, we cal- 
culate the time dependence of M shown in Fig|3l The 
oscillations of X{t) results in small x and y components 
of the magnetization. The magnetic moment precesses 
with a velocity Vx{t) = dX{t)/dt in the presence of the 
SSE-generated spin current {vx{Q) 7^ 0). We note that 
the maximum deflection angle X^ax depends implicitly 
on the spin current dynamics through the parameter a, 
as determined by eq. ([5]). Also the Fermi energy enters 
through the k dependence of a. The dynamics is a mix- 
ture of anisotropy-dominated precession and damping. 

Displacement- Let us consider the initial localized mag- 
netic moment distribution 

n — [to|| {t) sin C, ray{t), m|| {t) cos Q 

( = cos"^ [t£inh^{x/w)] (17) 

where w stands for the extension of the localized moment 
and ™||(0) — 1, my{0) = 0, a;(0) — 0. As concluded from 
Fig HI the moment is set in motion when subjected to 
the spin current. The velocity changes from positive to 
negative, which is different from the motion of a single 
Neel walP^ (the velocity decreases to zero in a fraction 



of a nanosecond, and the DWs stops completely.). 

Remarks and conclusions.- Our main result is that the 
SSE generated spin current in a wire may scatter from a 
localized magnetic structure setting it in an precessional 
and a displacement motion. The scattering also leads 
to a redistribution of the spin current in the wire and 
hence changes the ISHE signal. From these results con- 
clusion can be made on the influence of a collection of 
non-interacting localized moments but no statement can 
be made when they interact or even form clusters. We 
also note that the present conclusions do not apply to a 
domain wall (DW), (except for very close non-resonant 
(transversal) 180° DW pair, e.g. as in2I). In fact, our ini- 
tial finding^^ is that a single sharp DW is less affected by 
the spin current because the spin-current torques acting 
from left and right of the DW partially compensate. This 
is not so for an adiabatic or asymmetric DW because the 
T-gradient modifies the DW along VT. As for the exper- 
imental observation of the magnetic moment dynamics, 
it should be noted that the temperature gradient has to 
be sustained on the time-scale of this dynamics; a fast 
(e.g., femtosecond) strong heat pulse is inappropriate for 
our (constant VT, linear response) study and may cause 
locally a longitudinal dynamics of M. 
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